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Abstract 



The density fluctuations of nuclear matter are studied within a mean-field 
model in which fluctuations are generated by an external stochastic field. The 
constraints imposed on the random force by the fluctuation-dissipation the- 
orem are analyzed. It is shown that in the proximity of the borders of the 
spinodal region the assumption of a withe-noise stochastic field can be reliably 
used. The domain distribution of the liquid phase in the spinodal decompo- 
sition of nuclear matter is derived. The related distribution of fragment sizes 
compares favourably with the experimental fragment distribution observed in 
heavy ion collisions. 
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I. INTRODUCTION 



Semiclassical kinetic equations for the one-body phase-space density provide a powerful 
tool for studying the dynamics of complex processes occurring in heavy ion collisions [|] |3|] . 
However, these equations in their original version give a deterministic description for the 
evolution of the one-body phase-space density and their solution represents the mean value 
of this density at each time. Thus, they are not able to account for phenomena such as 
the nuclear multifragmentation observed in heavy ion collisions. In this process fluctuations 
about the mean phase-space density are believed to play an essential role (for a review 
on nuclear multifragmentaion see, e.g., Ref. ||). In the last decade, to remedy for this 
drawback, an extension of the transport theory has been proposed • This approach, that 
incorporates a fluctuating stochastic term into the kinetic equation, is usually known as the 
Boltzmann-Langevin equation, and was originally applied to the treatment of hydrodynamic 
fluctuations in the theory of classical fluids ||. In Refs. |5|-|7j collisions between nucleons 
in the nuclear medium are regarded as random processes and the diffusion coefficient of the 
Langevin (fluctuating) term is ultimately related to the amplitude of the nucleon-nucleon 
scattering. The last step is a particular case of the fluctuation-dissipation theorem. More 
recently, a new method to take into account fluctuations has also been proposed ||. In 
the approach of Ref. || the statistical fluctuations of the one-body phase-space density are 
directly introduced by assuming local thermodynamic equilibrium. The white-noise nature 
of the stochastic term is the basic assumption, generally shared by all works on this subject. 



The authors of Ref. [10], instead, have introduced an extension of the Boltzmann-Langevin 



theory by including a coloured-noise term in the stochastic force, the occurrence of such a 
term has been ascribed to the finite nucleon-nucleon collision time. Actual applications to 
nuclear problems of this interesting approach still have to be made. 

In the present paper, we study the density fluctuations and their time evolution by in- 
troducing a self-consistent stochastic field acting on the constituents of the system. The 
self-consistency condition is provided by the fluctuation-dissipation theorem. The evolu- 
tion of the fluctuations is treated within a linear approximation for the stochastic field. 
For simplicity we consider an infinite homogeneous system. First, without introducing any 
particular assumption, we prove that a withe-noise stochastic field cannot satisfy the self- 
consistency condition in general. Then, with reference to infinite nuclear matter, and within 
a collisionless mean-field approximation, we specify the particular conditions under which 
the withe-noise assumption for the stochastic field can be retained. These conditions are 
fulfilled for values of density and temperature lying in the proximity of the boundary of the 
spinodal region in the phase diagram, both inside and outside this region. Thus we consider 
nuclear matter in this physical situation and are able to solve the stochastic equation for the 
density fluctuations in a closed form. 
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With respect to the previous works on this subject here we consider a different 

source of fluctuations: Landau damping. This source is present even when collisions have a 
negligible role in the evolution of the system. 

In Sec. II we propose a procedure to determine the structure of domains formed within 
the system during a spinodal decomposition. If the fragmentation phenomenon observed in 
heavy ion collisions can be ascribed to a spinodal decomposition of the bulk of nuclei |TI , we 
are allowed to identify the pattern according to which nuclear matter is decomposing, with 
the fragmentation pattern, and can compare the results of our nuclear matter calculations 
with the fragment distribution observed in heavy ion collisions. This comparison is made in 
Sec. III. Finally in Sec. IV a brief summary and conclusions are given. 

Many papers, both theoretical and experimental, have been devoted to the multifrag- 
mentation problem. Here we mention only a few theoretical ones representing different 
approaches. In the statistical models of Refs. JT2.13], a complete statistical equilibrium of 



all degrees of freedom is assumed in a freeze-out volume and the various exit channels are 



sorted according to their statistical weight in the microcanonical ensemble. In Ref. |14| in- 
stead nuclear multifragmentation has been described in terms of "reducibility" and "thermal 
scaling" . This means that fragments are emitted practically independently of each other and 
the one-fragment probability is given by a Boltzmann factor. In the dynamical models of 
Refs. |5H7| clusters are constructed from the one-body phase-space density governed by the 



Boltzmann-Langevin equation [|15|-|17|]. In order to take into account the quantal nature of 



the system and the requirement of antisymmetrization, the Quantum Molecular Dynamics 
model |18| and the more sophisticated Fermionic (Antisymmetrized) Molecular Dynamics 



model P^J2D|] have been developed. In addition, percolation ET| and lattice-gas models 



22,23[ have also been introduced. These models are particularly suitable to deal with the 



critical phenomena which can be expected to occur in multifragmantation. To conclude 



this non-exhaustive survey, we mention the calculations of Ref. |24| that are based on the 
classical nucleation theory. 

Even if it should eventually turn out that multifragmentation must be ascribed to very 
complicated processes, we think that in any case our present approach could give some 
insight into the underlying mechanism. 



II. FORMALISM 

A. White noise assumption 

The mean-field approximation allows us to obtain a self-consistent equation for the time 
evolution of the one-body density. We assume that the time scale of the terms neglected in 
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the mean-field approximation is shorter than the characteristic times of mean-field dynamics. 
In order to take into account thermodynamic fluctuations, quantum effects and short-range 
correlations, we add to the mean-field a stochastic term similar to the random force in the 
Langevin equation. We assume that this additional field is a gaussian white noise with 
vanishing mean. In this case the time-evolution of the density is a markovian process. 
We will determine the conditions in which the white-noise assumption can be considered 
reasonable. 

The additional stochastic mean field will induce density fluctuations with respect to the 
mean density. To be more specific, we assume that at the time t = in the system is 
present a density fluctuation 8g(r,t = 0), with Sg(r,t) = g(r,t) — g , ( g is the density of 
the reference homogeneous state i.e the state towards which the system relaxes ). Within 
a linear approximation for the stochastic mean-field the Fourier coefficients of Sg(r, t) for 



t > are given by ( see for example Ref. p5| , Sec. 15 I ) 



Sg k (t) = Sg k (t = 0) - *ff ° } t D k (t - t') dt' + f D k (t - t')B k (t')dW k (t') , (2.1) 

JJk{uJ = U) Jo Jo 

where D k (t — t') is the response function of the nuclear medium and Dk(ui) its time Fourier 
transform. For symmetry reasons D^it — t') and its Fourier transform depend only on the 
magnitude of the wawe vector. In the second integral B^t^dW^t') gives the contribution of 
the stochastic field in the interval dt 1 . The real and imaginary parts of the Fourier coefficients 
IVk(i') are indipendent components of a multivariate Wiener process ||26|| . The fact that the 
stochastic field is real requires B£.(t) = B_ k (t) and W£(t) = W^(t). 

The stochastic part of the mean field is completely determined once the coefficients B^it) 
are known. In order to gain information about these coefficients we concentrate on the cor- 
relations of density fluctuations at equilibrium. Due to the independence of the components 
of the multivariate Wiener process Wy_(t), only the terms with k' = — k survive. Within a 
linear approximation, these correlations can be expressed by means of the same quantities 
that appear in Eq. (|2.1|) . This does, in a sense, correspond to the Onsanger hypothesis about 
the decay of deviations from equilibrium . The equation for the equilibrium fluctuations 



is obtained from Eq. fl2.1|) by moving the initial time to — oo, without including any particu- 
lar condition at finite times. Then the correlations are given by the equation ( the brackets 
denote ensemble averaging ) 

dhDkit-t^Dkit' - ti) S k (*i)S_ k (*i) . (2.2) 

-oo 



The time-translation invariance of the left-hand side of Eq. ( |2.2|) requires that the coef- 
ficients B]s_(t) must be constant. Taking the Fourier transform, Eq. ( |2.2j ) gives 
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= 2it5(uj + uj') D k (uj)D k (uj') |fik| 2 • (2.3) 

By exploiting the fluctuation-dissipation theorem 

2 

< (<^k^-k)M >= - _ e _pJ mD k {u) , (2.4) 

where /3 — 1/T is the inverse temperature, ( we use units such that h = c = ks — 1 ), we 
obtain for the coefficients the equation 

■ o|2_ 2 ,0^ 



We have used the relation D k {—uj) = D* k {u). Equation ( |2.5[ ) can be satisfied only if the 
right-hand side does not depend on uj. This can occur only in particular situations, thus the 
original white-noise assumption about the stochastic mean-field is not correct in general. 
This result is quite general, so we can conclude that for a perturbed system approaching an 
equilibrium state, fluctuations about the average trajectory cannot usually be accounted for 
by a white-noise stochastic force. 

Now, with reference to symmetric nuclear matter, we discuss particular physical situa- 
tions in which Eq.( [2.5|) can have a solution. Only in such conditions the assumption of a 
withe-noise stochastic field is valid. The relevant quantity is the linear-response function 
D k (uj). We evaluate D k (u) within a self-consistent mean-field approximation. In order 
to derive compact analytical expressions, here we use the linearized Vlasov equation for 
calculating the response function. This equation can be regarded as a semiclassical approx- 
imation to the random phase approximation, valid in the longwavelength limit. We also 
use a Skyrme-like form of the nucleon-nucleon effective interaction. Our self-consistent 
mean-field potential is given by 

U= a-^ + b ( _ d V 2 g , (2.6) 

Qeq Qeq 

where g eq is the saturation density of nuclear matter. For the parameters in Eq.( [2.6|) we take 
the values: 

a = -356.8 MeV, b = 303.9 MeV, a = - , 

6 



d = 130 MeV • fm 5 . 
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The values of a, b and a reproduce the binding energy ( 15.75 MeV ) of nuclear matter at 
saturation (g eq = 0.16 fin -3 ) and give an incompressibility modulus of 201 MeV. For the 
values of d we follow the prescriptions of Ref. . 



The response function is given by 



g »M = ; 7 nil : . ( 2 - 7 ) 

1 - AkDl >{ui) 

where Df, (ui) is the non-interacting particle-hole propagator, and 

A k = a— + -^q% + dk 2 (2.8) 

g eq a + 1 g^ 1 

are the Fourier coefficients of the effective interaction. Here go is the density of the reference 
homogeneous state. 

By substituting the expression ( [2.7D for Dkiyj) into Eq. (|2.5|) , we obtain 

2 _ 2 ImgfM 

This equation shows that the coefficients do not explicitely depend on the nucleon- 
nucleon effective interaction. However, we remark that D^\uo) is the propagator of inde- 
pendent particles that are moving in the mean-field of the reference homogeneous state, 

f/o= a— + M— ) Q+1 , (2-10) 

Qeq Qeq 

thus the interaction between constituents does enter, although not explicitly, into the ex- 
pression of £?k- 

We shall now show that in the classical limit uj/T <ti 1, the right-hand side of Eq.( |2.9| ) 
does not depend on ui, thus the assumption of a white-noise stochastic field can be considered 
valid in that limit. 

In the actual physical situations considered in this paper the values of the temperature 
are small enough with respect to the Fermi temperature so that the Pauli principle is still 
operating. Therefore the strength of the particle-hole excitations having energies much 
higher than kvp ( vp is the Fermi velocity ) can be considered negligible. Moreover the 
relevant values of the wave vector k turn out to be such that the quantity kvp is of the same 
order of magnitude as T. Thus the limit uj/T 1 also implies lojkvp 1. 



6 



The non interacting particle-hole propagator (u) acquires a very simple form in the 
longwavelength (Vlasov) limit. The imaginary part is 

^(.)=lm^ fde r ^e(l-^-), (2.11) 

where 



An k J de p V kv J 



h p g/3(6 p -A) _|_ i 



is the mean occupation number of nucleons with kinetic energy e p = p 2 /2m, and v = p/m. 
The effective chemical potential ft is measured with respect to the uniform mean field Uq. 

For ujjkvp <C 1 the imaginary part of D^(u>) is given by 

while the real part of Df\uj) in the longwavelength limit takes the form 

ReDf>(.) = -- 2 J *^ ( - 1 + - - .n — ^ ) . (2.13) 

For T sufficiently low with respect to ft, the most important contribution to the integral in 
Eq. fl2.13|) comes from a small domain of e p around ft. So we can take u/kv <C 1 in evaluating 
the integral, and obtain 

ReDf\u)=^+0((-^r). (2.14) 



dfx V kvp 



With D£\u) given by Eqs.(|2TT2|) and ( gig ), the right-hand side of Eq.QZ^) is indepen- 
dent of to the lowest significant order in cu/T. Thus, for u/T <C 1 the magnitude of the 
coefficients 5k is given by 



,2 . 2 2/^2 T 1 



|£ k M = -m^^-V— 5= f. (2.15) 

1 1 7T V <W e-^ + lfc v ; 

The phases of B^, instead, remain unknown. However, we will see that only the quantities 
\B\f\ 2 are needed to determine the probability distribution of density fluctuations. Finally, 
we remark that \B^\ for a given k, is determined solely by the density and temperature of 
nuclear matter. 
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The white-noise assumption is justified if the excitation strength is concentrated in a 
narrow range of energy close to zero. This condition requires that lmD k (u) is a sharply 
peaked function in the proximity of u = and is negligile elsewhere. The imaginary part 
of D k (u>) displays this feature for values of temperature and density near the borders of 
the spinodal region, since the pole of D k {u) lying on the imaginary axis, moves towards 
uj = as the system approaches the mechanical instability. This is shown in Fig.l, where 
we report lm.D k (u) calculated with Eq.fl2.7|) using the complete expression of D { k \u). With 
our effective interaction, for T = 5MeV the spinodal region starts at g c = 0.617 g eq . The 
values of go used in Fig.l are close to this critical value. 



B. Distribution of fluctuations. 



We now derive from Eq.(27T) the probability distribution for Sg^yt) in the limit u/T <C 1, 
and for values of temperature and density in the proximity of the spinodal region. The 
response function D k {uj) has a pole in the lower part of the imaginary c<j-axis, at a position 
given by 

iT k =i^ 2 (l + e-^^- '-k. (2.16) 

m z op, Ah 



Ao, (2.17) 



We have used the relation 

9fi _ d 2 f 
dg t dgl t 

where / is the free-energy density and Ao = A k =o- 

d 2 f 

In Eq.( |2.16D the relevant quantity is the isothermal stiffness — -j \t, which vanishes on 

OQq 

the boundary of the spinodal region. Since we limit our calculations to the proximity of the 

d 2 f 

spinodal region, we neglect tt^It with respect to Ao in evaluating DJt). Furthermore, in 

actual calculations the typical values of k which come into play are such that the term d k 2 
is smaller than |^4o|, thus we also neglect this term with respect to |^4o|- This approximation 
is consistent with the longwavelength limit adopted in the calculation of Dk(ui). 

Substituting into Eq.( [2.1|) the response function Dk(t — t') calculated with these approx- 
imations, the equation for the fluctuations 5g^(t) becomes: 

5&t(t) = 5g k (t = 0)e r ^ + B k f e Tk ^ dW k (t') , (2.18) 

Jo 
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where 



\Bi 



/2ttT 



1 + e-M) k , 



(2.19) 



and Tfc is given by Eq. (|2.16|) , neglecting the term dk 2 in Ak- We recall that Tk is negative, 

so that \Tk\ represents the damping rate of fluctuations, that vanishes for long wavelengths 

i d 2 f, 
when — — — t 



dgl 



0. 



Equation ( 2.18 ) represents an Ornstein-Uhlenbeck process [2B| with \Tk\ as drift coef- 
ficient and Bk as diffusion coefficient. The corresponding Fokker-Planck equation for the 
probability distribution P[5g k (t)] reads 



d 



1 , ~ 



5g k (t)P[dg k (t)]+-\B, 



2 



dSgl(t) 



[2.20) 



For simplicity we assume the state of the system at t = to be homogeneous on average 



( < Sgk(t = 0) >= for k 7^ ). Equation fl2.18|) says that this property holds during time 
evolution. In this case the solution of Eq.( [2.20|) is a gaussian distribution with zero mean 
value. Whenever it is necessary, a non vanishing mean value can easily be introduced. The 
explicit expression of the distribution P[8gk(t)] is 



P[5g k (t)] = N x e 
with the variance cr|(i) given by 



2 k CT fcW 



(2.21) 



4(t) = 4(t = 0) e 2 ^ + 



T 



[l-e 1 



fl' + dk 2 

Here the constant iV\ is a normalization factor and we have introduced the abbreviation 

d 2 f 
f = Z-L\ T 



(2.22) 



For t — > oo Eq.( 



dg 2 Q 

reproduces the usual gaussian approximation with variance 

T 



a k (t = oo) 



f" + dk 2 



(2.23) 
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for the equilibrium thermodynamical fluctuations pqj . We furthermore remark that 
Eq. ([2.22|) for the time evolution of the variance is similar to that obtained with different 
approaches in previous works on this subject [f7|,[2"9|. 

For later purpose we report also the distribution of the fluctuations in ordinary space: 
P M r,t)] = , iC l/«'«'.')*(ry,<)W, t ) (2m) 

where 

M(r,r',t)= ^£^ r - r,) - 



and A^2 is an appropriate normalization factor. 

The diffusion coefficients of Eq. (|2.19|) are derived by means of the fluctuation-dissipation 



theorem, which concerns only fluctuations about equilibrium. In Ref. |3(J a way has been 
suggested to extend the treatment of stable cases to processes where instabilities can develop. 
Following that suggestion we include in our approach the case of nuclear matter merged in 
the spinodal region. In practice, we still assume the validity Eq. (|2.21|) for the probability 
distribution of the density fluctuations, with the variance cr^(t) of Eq. ( |2.22| ) calculated with 
the values of temperature and density of the new situation. This amounts to treating the 
diffusion coefficients for the unstable case as an analytic continuation of the stable-case 
coefficients in the (g,T) plane. The reliability of such a procedure lies in the fact that 
both the growing rate and the diffusion coefficient B k change smoothly when the system 
crosses the stability boundary and enters the spinodal region. The pole of D k {oj) in turn, 
continuously moves along the imaginary axis from the lower part to the upper part of the 
complex cu-plane ( see Eq. fl2.16|) ). In order to preserve causality, the integral for calculating 
the Fourier anti transform Dk{t) must be performed along a path which cuts the imaginary 
axis above the pole. 



In the unstable case, the time behaviour of the variance <j\(t) in Eq. (|2.22|) is similar to 



that predicted by linear theories of the spinodal decomposition of alloys and fluids ( for an 
extensive review on this subject see Ref. |3]J ). The variance grows exponentially for the 
fluctuations with wave number 



fc<A:c= v?' (2 ' 25) 

while it tends to the aymptotic value a k {t = oo) of Eq. (|2.23| ) for k > k c . In particular the 
growth rate presents a maximum for k = ku = k c j\f?>. This means that the pattern of 
the regions which contain coherently correlated fluctuations is asymptotically characterized 
by the wavelength Xm = ^I^m- These features for the growth rate of unstable modes are 
analogous to those obtained in Ref. [[TIJ within a different scheme. 
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C. Size of fragments 



Starting from the probability distribution for density fluctuations given by Eq. (|2.24|) , 



we can determine the corresponding distribution for the size of the correlation domains. It 
has already been recalled that the stable and unstable cases can be treated within the same 
scheme. Thus we shall invesigate two different situations which could be explored by nuclear 
matter during a nucleus-nucleus collision: in one case the system is in the metastable region 
and relaxes towards a local minimum of the free energy, while in the other case the system is 
merged in the spinodal zone and develops density fluctuations which grow with time and will 
eventually lead to decomposition. According to our approximations, we limit our analysis 
in both cases to values of temperature and density in the proximity of the borders of the 
spinodal region. Moreover, we consider homogenous nuclear matter in both cases, and still 
assume that < 5gk(t = 0) >= for k ^ 0. 

Before performing explicit calculations, we make a few remarks. It is known that linear 
theories are unable to describe the late stages of the spinodal decomposition of alloys and 



fluids ( see Ref. |3T| and references quoted therein ). In particular, they predict a limiting 
value for the length scale that characterizes the pattern of the correlation domains. This 
value is given by the wavelength Xm for which the growth rate of fluctuations has a maximum. 



Instead, Monte Carlo simulations and experimental results [32] show a continuous coarsening 



of the domains with increasing time. However it has been argued in Ref. |33| that the early- 
time Monte Carlo results are consistent with a linear theory, provided that a stochastic force 
is included. 

In the physical situations considered in the present paper (heavy ion collisions), the value 
of the characteristic wavelength Xm is larger than 10 fm, beyond the size of the nuclear system 
involved. Moreover, the corresponding growth time ^-/Tk M is of the same order of magnitude 
as the characteristic times of the nucleus-nucleus collisions in the energy range considered 
here. Thus the fluctuations with wave number ku are still far from being the predominant 
ones in this time interval. This means that the processes that we are investigating correspond 
to an early stage of the spinodal decomposition. Then we can expect reliable results from 
our approach, at least at a qualitative level. 

From Eq.( |2.24| ) we obtain the usual expression for the equilibrium correlation function 

G{\v - r'|) = — - e - -, (2.26) 

An a |r — r | 

where 
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is the correlation length. This quantity, which represents the average extension of the cor- 
relation domains, can be obtained by an appropriately weighted integral of the correlation 
function: 

f= J drdr'F{r,r')G(\r-r'\). (2.27) 

The function F(r, r') is a suitable weight function. Here we extend this relation between 
averaged quantities to fluctuating quantities, for systems both at equilibrium and out of 
equilibrium. We then assume that the size of correlation domains at time t is given by a 
quadratic functional of the fluctuations Sg(r,t): 

H) fdrdr'F{r,r')G{\r-r'\,t) ' 1 ' 

where L(t) = < b > is the length scale that characterizes the pattern of the domains, and 
G(\r — r'\,t) is the correlation function for systems out of equilibrium. The latter quantity 
is the space Fourier transform of the variance of Eq. (|2.22[) . 

In order to simplify calculations, we further choose for F(r, r') a separable form. The 
requirement that b should be positive for any function 5g(r,t), enforces a symmetric form 

F(r,r')= /(r)J(r'). (2.29) 

of the weight function. This form allows us to obtain a closed expression for the probability 
distribution of b. In addition, with this choice the final results are entirely independent of 
the function /(r). 

Now we derive the probability distribution for b as a function of the length scale L(t). 
Later we shall give a procedure for determining L(t). 

For a given probability distribution P[8g(r,t)], the related probability distribution for b, 
at a given time t, can be obtained by means of the functional integral 

P(M) = / d[Sg(r,t)]6(b- ^ J dr5g(r, t)F(r)5g(0, i?)P[6 e (T, t)} , (2.30) 
where we have put 

C= J dvdv'f{v)G{\v-v'lt)f{v') (2.31) 

in order to simplify the notation. 

With the distribution of Eq. (|2.24 ) and using the integral representation of the 5-function 
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5(x) = / ^e^ x 

2lX J 



the equation for the distribution P(b, t) takes the form 
N> 



P(M) = [ drje^ 

(2.32) 

-i | dr dr'^(r, t) (M(r, r', t) + 2irj^-f(r)f(r')) 5q(t', t) 



x / d[8g(r,t)]e 



The functional integral is of gaussian type and allows us to express the result of the integra- 
tion in closed form: 



P(b,t) = — / d^'l" — ; TTTT^ — r • (2.33) 

The quantities M and F are infinite-dimensional operators, with matrix elements M(r, r', t) 
and /(r)/(r') respectively, in the coordinate representation. 

The determinant in the last equation can be factorized as 

det^- [M + 2irj^-F] = detAdet[l + 2^^M _1 F] , (2.34) 
27r C 2n C 

where 1 is the unit matrix. The square root of the first factor on the right-hand side and the 
normalization constant N2 of Eq. ( |2.33| ) coincide and cancel. What remains to be evaluated 
is the inverse of the square root of the second determinant. For this purpose we write the 
determinant in exponential form and expand the exponent in a power series. Thus, we obtain 
the following formal expression 

L(t) . 1 --Trln(l + 2in^-M- 1 F) 

(det[l + 2t V =^-M- 1 F])-2 = e 2 1 1 C ' 

(2.35) 

ly I(_i)i+«( 2 ^^)™Tr(M~ 1 F) n 
2 n C ' 



e 



We recall that the matrix element M -1 (r, r', t) and the correlation function G(\r — r'|,£) 
coincide. Because of the separable form chosen for the function F(r, r'), Eq. (|2.29|) , the trace 
operation on the generic n-term of Eq. fl2.35|) simply yields C n . Thus the series can be 
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resummed and gives ln(l + 2irjL(t)). 
form 



Then the probability distribution P(b, t) acquires the 



P(b,t) = — [d V 



e ir]b 



(1 + 2ir ] L(t)) 1 2 



A simple integration in the complex rj plane gives the final result 



P(b,t) 



1 1 



e -b/(2L(t)) 



(2.36) 



^2L(t)b 



From the probability distribution of the domain size we can derive the distribution of the 
number of nucleons A that are contained in a correlation domain, assumed to be spherical. 
For a homogeneous liquid the relation between A and the size b is b = 2roA 1 ^ 3 , where r is 
determined by the actual density. With a simple transformation of variables we obtain for 
the probability distribution of A, P(A,t), the equation 



Further, to take into account that A is a discrete variable we express the probability of 
finding a correlation domain containing A nucleons, Y(A), through the integral 



The distribution P(A,t) and the probability Y(A) are completely determined once the 
ratio between the length scale L(t) and the mean interparticle spacing r$ is fixed. The pa- 
rameter L(t) sets the scale for the decrease of the correlation function G(r, t) with increasing 
r. We can obtain an estimate of L(t) by analyzing the behaviour of G(r,t) as a function of 
r at a given t. The correlation function is initially determined by the variance a\{t = 0), 
then, in the stable case, it asymptotically assumes the form given in Eq.( [2.26j ), with the ap- 
propriate correlation length £ = L(t = oo), while in the unstable case, it acquires a damped 
oscillatory behaviour carachterized by the asymptotic wavelength X M . In order to illustrate 




(2.37) 




(2.38) 



For large A, Y{A) tends to coincide with P(A,t). 



III. RESULTS 
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the general features of the function L(t), we simply assume that the initial fluctuations are 
negligible, o\{t = 0) ~ 0. In this case the function G(r,t) is completely determined by the 
density and temperature of nuclear matter. Here we consider two sample values for the den- 
sity ( f?o — 0.65^ eg and g = 0.58^ e? ) and a single value for the temperature ( T = 5 MeV ). 
This temperature is in the range of values expected for the nuclear multifragmentation pro- 
cess [|J. The two corresponding points in the phase diagram (g, T) lie in the metastable 
region and in the spinodal region respectively, and are sufficiently close to the boundary of 
the spinodal zone to justify our assumption of a white-noise stochastic field. In Figs. 2 and 
3 we show the behaviour of G(r,t) as a function of r at three different values of time, both 
in the stable and unstable situations. In the stable case of Fig. 2, a simple inspection of the 
behaviour of G(r, t) shows that it is reasonably well reproduced by a function like that on the 
right-hand side of Eq.( |2.26| ) (obviously with £ replaced by L(t) ). We adopt such a form for 
G(r,t), then, by comparison with its true behaviour shown in Fig. 2, we can determine L(t). 
For the unstable case shown in Fig. 3, the situation is slightly more involved because the 
asymptotic regime is reached only after a very long time. For this case, we simply assume 
that L(t) does coincide with the distance at which the value of G(r,t) is reduced by 80% 
with respect to its value at r = 1 fm (because of our approximations we cannot expect the 
present approach to be reliable for distances shorter than lfm). 

At a given time t the value of the length scale L(t) depends strongly on the distance from 
the boundary of the spinodal zone, the shorter this distance, the larger is L(t). In Fig. 4 the 
calculated length L(t) is displayed as a function of t for the two chosen sets of parameters. 
The values of t are in the range that is relevant for nuclear fragmentation Figure 4 
shows that for t ~ 200fm/c, L(t) pratically reaches its asymptotic value (L(oo) ~ 3.0 fm) 
in the metastable situation, whereas in the unstable case L(t) is still much smaller than 
L(oo) w 12 fm. 

In the two physical situations considered here, two different processes could drive nuclear 
matter towards a spinodal decomposition. In the metastable case, if the density fluctuations 
are large enough, the nuclear system can explore the unstable region for a time sufficiently 
long to move towards a phase separation. In the unstable case instead, fluctuations grow 
with time until they cause the decomposition of the nuclear system. In both cases we expect 
that the pattern of domains containing the liquid phase is determined by the probability 
distribution P(b,t) or P(A,t) of Eqs.flESSD and ( g37fl . 

In order to assess the degree of validity of our approach, we compare the results of our 
calculations with the corresponding experimental data by identifying the probability Y(A) 
of Eq.( |2.38|) with the distribution of the fragment yield. Since experimentally the fragments 



are detected according to their charge, we have to transform P(A,t) and Y(A) into the 
corresponding functions of Z. We assume a homogeneous distribution also for the charge 

Z = - -A, with a = (N — Z)/A, and use a = 0.2, which corresponds to the average 
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asymmetry of the nuclear systems considered. 

In Fig. 5 the probability Y(Z) is displayed as a function of Z on a double logarithmic 
scale for three different values of the ratio L(t)/rQ. The range of values for L(t)/r has been 
chosen in accordance with that of L(t) in Fig. 4. Figure 5 shows that Y(Z) can be fit with 
good accuracy by a power law Y(Z) = Y Z~ Te ff . The values of the effective exponent, r e ff, 
lie between 1.17 for L(t)/r = 4 and 1.42 for L(t)/r = 2. 

The power-law behaviour of the fragment yield and the determination of the exponent 
have been the subject of several experimental studies of multifragmentation ( see for example 
the recent papers [|34l , j35|l ). The observed values of the exponent are in the interval ~ 1.2 — 1.5 
for nuclear rections with beam energies lower than ~ 40 A MeV, whereas they exceed the 
value of 2 at higher energies |34]j35|1 . A value of the exponent r e ff > 2 can be unlikely 



reproduced by our calculations because we would need an unreasonably low value for the 
ratio L(t)/r . However, in various papers [35-37] it has been remarked that the effects of 
collective motions, that have not been taken into account by our present approach, should 
become more important with increasing beam-energy. 

Figures 6 and 7 show a comparison between the charge distributions predicted by our ap- 
proach, Y(Z), and recent experimental data obtained by the Multics/Miniball collaboration 
for Au + Au collisions at an incident energy of E = 35 A MeV |38j and by the INDRA Col- 
laboration for 129 Xe + Sn and 155 Gd + 238 U collisions at E = 32 A MeV and E = 36 A MeV 
respectively []39 |. The calculations have been performed for three values of the parameter 



L(t)/rQ. We have normalized the experimental distributions to one in order to perform the 
comparison on an absolute scale. We can see that the agreement between experimental data 
and the calculated charge distributions is quite satisfactory for Z < 30 4- 35 and that for the 
lighter fragments the experimental points are better reproduced with larger values of the 
ratio L(t)/r Q . For Z > 30 4 35 the observed distribution presents a slope steeper than that 
predicted by our calculations. This faster decrease should be ascribed to finite-size effects 
which have not been included in our nuclear matter treatment. 



IV. SUMMARY AND CONCLUSIONS 



We have studied the density fluctuations associated with a one-body treatment of nuclear 
dynamics. In our approach the fluctuations are generated by adding a stochastic term to 
the mean field. This additional random force is determined by a self-consistency condition 
required by the fluctuation-dissipation theorem. We have treated the effects of the stochastic 
field in linear approximation and this has allowed us to express the time evolution of the 
fluctuations in a closed form. 

First we have analyzed the nature of the stochastic field and have shown that in general 
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a white-noise assumption for the stochastic field is not consistent with the fluctuation- 
dissipation theorem. Then we have studied the particular physical conditions in which the 
white-noise nature of the stochastic term can be retained. These conditions include hot 
nuclear matter at a temperature T 5 MeV, where the system can be still considered 
degenerate. We have found that for a Fermi system the treatment of density fluctuations 
by means of a white-noise stochastic term is justified when the limit uj/T <C 1 gives a 
reasonable approximation to the density-density response. This condition is better satisfied 
when the density and temperature of the system are close to the borders of the spinodal 
region in the (g,T) plane. Thus, in the limit uj/T <C 1 the equilibrium fluctuations can be 
adequately described by means of thermodynamic functions and we can expect that in this 
limit the purely quantum fluctuations will play a negligible role also for systems not too far 
from equilibrium. We have extended the results obtained for the probability distribution 
of a metastable system to unstable situations. This has been achieved by extrapolating 
the relevant quantities across the boundary of the spinodal region. Because of the linear 
approximation used for evaluating the response of the system to the stochastic force, the 
fluctuations have a gaussian probability distribution. 

In the final part of this paper we have introduced a procedure to determine the size 
and mass distributions of the domains containing correlated density fluctuations, then we 
have compared the obtained mass distribution to the yield of light fragments observed in 
the multifragmentation of heavy nuclei. The procedure proposed here is quite general and 
can be applied to any gaussian fluctuation distribution 

Our approach can account both for the observed power-law distribution and for the 
value of the effective exponent found experimentally, but for the exponent the agreement is 
limited to collisions with beam energies lower than ~ 40 A MeV. This discrepancy between 
our predictions and the observed values of the effective exponent in collisions of higher 
energies deserves further investigations. A detailed comparison with experiment has shown 
that our approach fairly reproduces the measured charge distributions for Z < 30 4- 35. 
Since we are dealing with infinite nuclear matter, we expect to overestimate the number of 
fragments having a large fraction of the mass of the emitting source. 

Finally, we remark that the obtained mass distribution contains only one parameter, the 
ratio between the time-dependent length scale of domains L(t) and the mean interparticle 
spacing r . This ratio can become large. A more detailed comparison of the present model 
with experimental data could also give an estimate of the time required by the system to 
break up. 



17 



REFERENCES 



[1] G.F. Bertsch and S. Das Gupta, Phys. Rep. 160, 189 (1988). 

[2] P. Schuck, R.W. Hasse, J. Jaenicke, C. Gregoire, B. Remaud, F. Sebille, and E. Suraud, 

Prog. Part. Nucl. Phys. 22, 181 (1989). 
[3] A. Bonasera, F. Gulminelli, and J. Molitoris, Phys. Rep. 243, 1 (1994). 
[4] B. Tamain and D. Durand, in Proceedings of the International School Trends in Nuclear 

Physics, 100 years later, Les Houches, Session LXVI, edited by H. Nifenecker, J. -P. 

Blaizot, G.F. Bertsch, W. Weise, and F. David (North-Holland, Amsterdam, 1998), p. 

295. 

[5] S. Ayik and C. Gregoire, Nucl. Phys. A513, 187 (1990). 

[6] J. Randrup and B. Remaud, Nucl. Phys A514, 339 (1990). 

[7] M. Colonna, Ph. Chomaz, and J. Randrup, Nucl. Phys. A567, 637 (1994). 

[8] M. Bixon and R. Zwanzig, Phys. Rev 187, 267 (1969). 

[9] M. Colonna, M. Di Toro, A. Guarnera, S. Maccarone, M. Zielinska-Pfabe, and H.H. 
Wolter, Nucl. Phys. A642 449 (1998). 
[10] S. Ayik and J. Randrup, Phys. Rev. C 50, 2947 (1994). 
[11] G. F. Bertsch and P.J. Siemens, Phys. Lett. B 126, 9 (1983). 

[12] D.H.E. Gross, Rep. Prog. Phys. 53, 605 (1990); D.H.E. Gross, Phys. Rep. 279, 119 
(1997). 

[13] J. P. Bondorf, A.S. Botvina, A.S. Iljionov, I.N. Mishustin and K. Sneppen, Phys. Rep. 
257, 133 (1995). 

[14] L.G. Moretto, R. Ghetti, L. Phair, K. Tso, and G.J. Wozniak, Phys. Rep. 287, 249 
(1997). 

[15] A. Guarnera, M. Colonna, and Ph. Chomaz, Phys. Lett. B 373, 267 (1996). 
[16] M. Colonna, Ph. Chomaz, and A. Guarnera, Nucl. Phys. A613, 165 (1997). 
[17] M. Colonna, Nucl. Phys. A630, 136c (1998). 

[18] J. Aichelin, Phys. Rep. 202, 233 (1991); R. Nebauer, J. Aichelin, and the INDRA 

Collaboration, Nucl. Phys. A658, 67 (1999). 
[19] H. Feldmeier and J. Schnack, Prog. Part. Nucl. Phys. 39, 393 (1997). 
[20] A. Ono and H. Horiuchi, Phys. Rev. C 53, 2958 (1996). 

[21] W. Bauer and A. Botvina, Phys. Rev. C 52, R1760 (1995); T. Gharib, W. Bauer, and 

S. Pratt, Phys. Lett. B 444, 231 (1998). 
[22] X. Campi and H. Krivine, Nucl. Phys. A620, 46 (1997). 

[23] S. Das Gupta, J. Pan, I. Kvasnikova, and C. Gale, Nucl. Phys. A621, 897 (1997). 
[24] L.V. Bravina and E.E. Zabrodin, Phys. Rev. C 54, R464 (1996). 

[25] L. E. Reichl, A Modern Course in Statistical Physics (Edward Arnold Publishers, Lon- 
don, 1980). 

[26] C.W. Gardiner, Handbook of Stochastic Methods (Springer- Verlag, Berlin, 1997), Ch. 4. 



18 



[27] W.D. Myers and W.J. Swiatecki, Nucl. Phys. 81, 1 (1966). 

[28] L.D. Landau and E.M. Lifshitz, Statistical Phisics 3rd Edition Part 1 (Pergamon Press, 

Oxford, 1982) Ch. XII. 
[29] D. Kiderlen, Nucl. Phys. A589, 320 (1995). 

[30] D. Kiderlen and H. Hofmann, Phys. Lett. B 353, 417 (1995); H. Hofmann and D. 

Kiderlen, Int. J. Mod. Phys. E7, 243 (1998). 
[31] J.D. Gunton, M. San Miguel, and Paramdeep S. Shani, in Phase Transitions and Critical 

Phenomena, edited by C. Domb and J.L. Lebowitz (Academic Press, London, 1983), 

Vol. 8, p. 267. 

[32] J.L. Lebowitz, J. Marro, and M.H. Kalos, Acta Metall. 30, 297 (1982). 

[33] J. Marro and J.L. Valles, Phys. Lett. A 95, 443 (1983). 

[34] T. Li et al, Phys. Rev. C 49, 1630 (1994). 

[35] M. D'Agostino et al, Phys. Rev. Lett. 75, 4373 (1995). 

[36] W.C. Hsi et al, Phys. Rev. Lett. 73, 3367 (1995). 

[37] G.J. Kunde et al, Phys. Rev. Lett. 74, 38 (1995). 

[38] P. Desesquelles et al, Multics/Miniball Collaboration, Nucl. Phys. A633, 547 (1998). 
[39] M.F. Rivet et al, Phys. Lett. B 430, 217 (1998). 
[40] M. D'Agostino et al, Phys. Lett. B 371, 175 (1996). 



19 



Figure captions: 



Fig.l Imaginary part of the response function D^.{u) of Eq. (|2.7| ) for hot nuclear matter 
(T = 5 MeV) at different densities approaching the critical value g c = 0.617 g eq (dotted 
line: g = 0.70 g eq , dashed line: g = 0.65 g eg , full line: g = 0.63 g eq ). The value of k 
is 0.1 k F . 



Fig.2 Spatial behaviour of time-dependent correlation function G(r,t) for nuclear matter 
on the stable side of the spinodal curve (go = 0.65 g eq , T = 5 MeV). The three curves 
correspond to different values of t (full line: t = 50fm/c, dashed line: t = lOOfm/c, 
dotted line: t = 200fm/c). 



Fig.3 Same as Fig.2, but for the unstable case (go = 0.58 g eq ). 



Fig.4 Behaviour of the length scale L(t) within the spinodal region (full line: g — 0.58 g eq ) 
and outside it (dashed line: go = 0.65 g eq ). 

Fig.5 Fragment distribution Y(Z) calculated for different values of the ratio L(t)/r Q . Full 
line: L(t)/r = 4, dashed line: L(t)/r = 3, dotted line: L(t)/r = 2. 

Fig.6 Comparison of fragment distribution Y(Z) calculated for L(t)/ro = 6 (full line), 4 
(dashed line), 2 (dotted line) with experimental distribution for the reaction Au + Au 
at 35 A MeV. Data from Ref. [081 have been normalized to one. 



Fig. 7 Same as Fig.6, but for the reactions 129 Xe + Sn at E = 32 A MeV (triangles) and 
155 Gd+ 238 U at E = 36 A MeV (circles). Data from Ref. [|39H have been normalized to 



one. 
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